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Abstract 

The physics of spontaneous chiral symmetry breaking in the case of the simultaneous 
presence of a magnetic field and a fermionic quartic interaction is discussed for both 
local and nonlocal kernels in 2+1 and 3+1 dimensions. The approach is based on the 
use of Valatin-Bogoliubov canonical transformations, which allow, in the absence of 
fermionic quartic terms, to completely diagonalize the Hamiltonian and construct 
the vacuum state. 

PACS: 12.38.Aw, ll.30.Rd, ll.lO.Kk 



1 Introduction 

During the past few years, the behavior of QCD and, more specifically, the 
problem of dynamical chiral symmetry breaking in the presence of a strong 
magnetic field have attracted some attention [1,2,3]. It has been shown that a 
constant magnetic field acts as a strong catalyst of dynamical chiral symmetry 
breaking, thus leading to the generation of a fermion mass [1,2]. The physics 
behind this effect is easy to understand: the motion of a charged particle is 
restricted to directions perpendicular to the magnetic field and this leads to 
a dimensional reduction {d ^ d — 2) in the dynamics of the fermion pair- 
ing. Different methods can be used to study this problem. Among them, the 
Schwinger proper time method [4] allows to obtain the exact expression for 
the fermion propagator S{x, y) in a constant magnetic field and evaluate the 
condensate in a limiting procedure, (0|^/'^/'|0) = — lim Tr S{x,y). In an alter- 
native approach, one can develop the propagator in terms of Landau poles [5] 
and obtain the condensate as the contribution coming from the lowest Landau 
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level (LLL). In this sense, the LLL plays a role similar to the Fermi surface 
in the Bardeen-Coopcr-Schrieffer (BCS) superconductivity theory. The draw- 
back of the above methods is the fact that they do not reveal explicitly the 
vacuum structure in terms of fermion (antifermion) creation and annihilation 
operators. On the other hand, not only calculations are simpler in the oper- 
ator formalism but also the interpretation of the pairing structure as well as 
its nonperturbativc nature are more transparent [6]. In this paper we adhere 
to the spirit of Ref. [6] and develop a general operator formalism necessary to 
study the vacuum structure and fermion condensation in systems unyielding 
to the usual mean-field approach. In particular, we study the mechanism of 
fermion condensation in models where, besides an external magnetic field, a 
strong local or nonlocal confining interaction is also present. 

The paper is organized as follows: in section 2 we introduce the formalism, 
based on the use of successive Valatin-Bogoliubov canonical transformations. 
In section 3 we exemplify its use by working out the complete diagonalization 
of the Hamiltonian of a free fermion in the presence of a constant magnetic 
field for both the (2+1)- and (3+l)-dimensional cases. Section 4 is devoted to 
test our formalism against the well-known results of the local Nambu-Jona- 
Lasinio (N JL) model [7] , setting the stage for the new results concerning non- 
local kernels in the presence of magnetic fields. In particular, we show how 
to translate the formulae and results obtained within the present formalism 
to the ones usually obtained through the use of mean-field techniques [1,2]. 
This translation is a validation of our formalism. The new and far more diffi- 
cult case of a nonlocal fermion kernel with a constant magnetic field is dealt 
in section 5. For simplicity, we have chosen a harmonic confinement. Within 
the present formalism other fermion kernels, notably the linear confinement 
kernel, could be considered as well. However, the linear kernel, which involves 
integro- differential mass-gap equations versus the simpler differential equa- 
tions proper of the harmonic confinement, puts, for the purpose of this paper, 
quite unessential calculation complexities. The end of this section is devoted 
to show that this harmonic confining, extended NJL model, already behaves, 
in the presence of a magnetic field, quite differently from the local-kernel NJL 
model. Our conclusions are presented in section 6. 



2 The approach 

The Lagrangian of a relativistic fermion in an external field ^4^ has the stan- 
dard form 
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where = + ieA^ is the covariant derivative. ^ Different gauges can be 
used to solve the Dirac equation in a constant external magnetic field. In 
what follows we choose the Landau gauge = —By , where B is the 
magnetic field strength. In 2+1 dimensions, the presence of a homogeneous 
magnetic field breaks polar symmetry explicitly, while in 3+1 dimensions, the 
spherical symmetry is broken. As is well known, this fact leads to the use 
of Landau levels as the appropriate basis for the remaining symmetry. The 
energy spectrum is given by E"^ = + pi + {2n + 1 — s)|e-B|, where s = ±1 
defines the spin orientation and is the momentum along the z-direction ^ . 
The discrete number n = 0, 1, 2, . . . describes the Landau levels. 

The Hamiltonian of the problem is given by 

Jd^xtlj\x){pm + a-p-ea-A)'ijj{x), (2) 

where D is the space dimension and we use the standard notation /3 = and 
a = '-f^'y. In order to construct the Dirac field in a magnetic field we start 
from the plane wave decomposition of the free field. 



^(^) = E ^7 {^M asp + v,ip) bl_A e^P-, (3) 



where V is the volume. Clearly, in 2+1 dimensions the spin quantum number 
s has no meaning and it should be omitted in all the corresponding formulae. 
The operators agp and bsp satisfy the usual anticommutation relations. The u 
and V spinors are the solutions of the Dirac equation for positive and negative 
energies, respectively. 

Since the quantum field of a particle in a magnetic field is usually developed 
in a basis indexed by the Landau levels n, a change of basis is needed in order 
to relate it with the plane wave solution. To achieve this, we first perform 
the canonical transformation which relates the spinors Us'ip) ,Vs'(p) with the 
momentum-independent (p — 0) spinors Ug ,Vs, 

-Eo:mH^\, (4) 

\Vs'(P) 

where Oss'ip) is the rotation matrix. The resulting field maintains the canon- 
ical structure (3). Next, in order to index the spinors and operators by the 





^ We recall that in 2+1 dimensions there are two inequivalent representations of 
the Dirac algebra, described by the matrices 7'* and —7'^, with 7° = 0-3 , 7^ = 
iai , 7^ = ZC72; Cj arc the Pauli matrices. The corresponding chiral version of 
the problem can then be formulated using the 4-dimensional spinor representation 
7^^ = diag (7'*, —7'*) . 

2 In 2+1 dimensions, Pz = and s = +1 (s = -1) ioi E > {E < 0). 
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Landau levels, wc expand the field in a complete set of Hermite polynomials 
Hn through the identity 



where 



OxPy 



27r^ eUJn{i)uJn 



yj ' 



n=0 



(5) 



(6) 



and 



(7) 



£ — \eB\ is the magnetic length. The corresponding change of basis is 
given by 



"snp 



^ fnss'iiPy) 



s',Py 



(Xg'p 



Us' 
V.I 



, (8) 



where p = {px,Pz) in 3+1 dimensions and p = Px m 2+1 dimensions, and 



the operators satisfy the usual anticommutation relations {a 



t „ 

snp 5 '^s' 



n'p' } 



{bsnp : bs'n'p'} = ^ss' ^nn' ^pp' ■ The matrices Tn{ipy) and Tn{^) can be obtained 
after a simple algebra. 

The last step consists in obtaining the mass gap equations which lead to the di- 
agonalization of the Hamiltonian. There are several approaches one can adopt 
to obtain the mass gap equation. It can be derived as the condition for the 
vacuum energy to be a minimum [8,9], or in the form of a Dyson equation 
for the fcrmion propagator [9], or as a Ward identity [10]. Here we shall de- 
rived it by imposing the vanishing of the anomalous (nondiagonal) terms in 
the Hamiltonian, corresponding to the direct creation (or annihilation) of a 
fermion-antifermion pair [11,12]. Aiming at this, we perform a second canon- 
ical transformation defined as 



"'snp 
y ^sn—p J 



I 



y^^Rnss'jPz) 



ds' 



\ 



np 



\ b\'n-p ) 




^nss' {P 




(9) 



where the rotation matrix Rnss'iPz) is obtained by requiring the anomalous 
terms in the Hamiltonian to vanish. 



After this diagonalization procedure, the vacuum state in the magnetic field 
can be easily constructed. Indeed, each Valatin-Bogoliubov transformation on 
the operators and spinors corresponds to a transformation of the vacuum such 
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that the latter is annihilated by the new operators. In the next sections we 
give explicit realizations of this construction. 

The approach presented above can also be extended to the study of systems 
having, besides an external magnetic field, a quartic fermion interaction of the 
form 

Hir^t^ J d^'x J d''yij{x)Tij{x)V{\x -y\)ij{y)Tij{y) , (10) 

where F denotes the vertex (Dirac, color) structure, and V is an attractive 
potential. This type of kernels has been extensively used in hadronic phe- 
nomenology [12,13]. In this case, we write the total Hamiltonian in normal- 
ordered form with respect to the new vacuum. Using the Wick contraction 
technique, the Hamiltonian can be decomposed in three terms, 

H ^ ■.Ho: + ■.H2: + -.H^:, (11) 



containing 0, 2 and 4 creation and/or annihilation operators, respectively. In 
the presence of a magnetic field, K — Pm+a-p—ea-A. The Greek indices refer 
to the two different Dirac fields, ipi and ^"2, corresponding to two inequivalent 
representations of the Dirac algebra in the (2-|-l)-dimensional case. In 3-1-1 
dimensions, such indexing is not necessary. 

The contribution to : : coming from H^nt can be diagrammatically expressed 
in 2-1-1 dimensions as 

Vlr^iV^^/'IrV'i + 4t^2v4^^2 4 [a -b]^, + 4 [c - ^2 , 
4r^|JlV42^i^2 + 42^i^2v4iTiJi 4 [d] ^1 + 4 [b] ^2 , (12) 
4iTiiJ2v42^^i + 4r^iV^Ir^2 4i [c] ^1 + 4 [a] ^2 , 

where A, B, C and D correspond to the diagrams presented in Fig. 1. On the 
other hand, in 3-1-1 dimensions the contribution is simply 

4T'iijv4^'il^ — ^4[a-b]'iIj. (13) 

Once the terms contributing to the quadratic part : H2 : have been identi- 
fied, its diagonalization proceeds through a sequence of Valatin-Bogoliubov 
transformations, as given in Eqs. (4), (8) and (9). 
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A C 



Fig. 1. Diagrams contributing to :i?2 : • The dashed hnes represent the interaction 
potential V and the single (double) lines correspond to the fermion propagator of 
the field ipi {ijj2) in the case of 2+1 dimensions. 

We note that, in the calculation of the diagrams of Fig. 1, the fermion propa- 
gators associated with ipi {i = 1,2) in 2+1 dimensions take the explicit form 



(14) 



while in the (3+l)-dimensional case the propagator is given by 

J{{x-x')px+{z-z')p^) 



l^^x^z n,p,s 



(15) 



3 Free fermions in a magnetic field 



As a first example of the applicability of our approach, we consider the problem 
of a fermion in 2+1 dimensions in the presence of a constant magnetic field. 
In this case, the rotation matrix 0{p) corresponding to the first canonical 
transformation (4) is simply given by 

I COS0 - sin {py + ip^)\ p \p\ 

0{p)=\ , P=r-|, tan20=^. 

ysin (f) {py -ip^) COS0 J \P\ 

(16) 

The vacuum associated to the operators dp and bp has the form 

|0) = l[(cos + sin alblp)\0) , (17) 
p 

so that OplO) = 0, bp\0) = 0. 
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The second canonical transformation defined in Eq. (9) is performed with the 
rotation 



Rn 



cos 6n — sin 6 
sin 9„ cos 9r, 



(18) 



where the angles 9n are determined by imposing the vanishing of the anoma- 
lous terms in the Hamiltonian. A simple algebraic computation yields the 
following mass gap equations, 



m sin 200 — , n — , 

m sin 2dn - i"^ V2n cos = , n > , 



which for any value of n have the solution 



v2n 

tan 29n = — j- 
mi 



(19) 



(20) 



It remains to construct the vacuum state in the magnetic field \0)b , which is 
annihilated by the operators Unp^ and bnp^ ■ We find 



1 0) B = n (cos + sin 9n alp^ bi_p^ ) 1 0) 



(21) 



n, Px 



At this point, it is worth commenting on the question of dynamical symmetry 
breaking in the presence of the magnetic field [2,6]. This problem is connected 
with the appearance of a fermion condensate b(0|^^|0)b 7^ in the limit 
where the mass parameter m goes to zero. Writing the field in terms of the 
new canonical operators one can easily obtain 

lim b(0|VSV'|0)b = - lim (cos 2^o + 2 V cos 2^^ ) 

= - lim — - 1 + 2 y , = , (22) 



i.e. the spontaneous breaking of the original U (2) symmetry of the Lagrangian 
(1) occurs even in the absence of any additional interaction between fermions [2] . 
This result is specific to the (2+1) -dimensional Dirac theory in an external 
magnetic field. 

Let us now consider the problem in 3+1 dimensions. The rotation matrix 
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Oss'ip) in Eq. (4) is given in this case by 



Oss'ip) 



^ cos Sss' — sin (j) Mss' ^ 
^smcf) Mlg, COS0 5ss' J 



m 



-(p-0-)(i(T2). 

(23) 



The new vacuum |0) can be generated from the |0) vacuum following the 
standard procedure and we find 

|6> = n f + ^ cl + i^cf ) 10). 



2 • 2 P 



(24) 



which obviously satisfies a^plO) = 0, 6sp|0) = 0. 

Next, we define the rotation matrix Rnss'iPz) of Eq. (9) as 



Rnss'iPz) — 



COS 6n Sss' -sinOnMnss' 



^ sin 6nMl^^, COS On 



^ — sin (fin COS (fin 

y COS cpn sincpn 



(25) 



with Mn chosen so that M^Mn = t and (po — 0. Once more, the angles 9n 
and (fin are determined by requiring the anomalous terms in the Hamiltonian 
to vanish. This leads to the following mass gap equations. 



m sin 200 — Pz cos 29q — , 



msin2^„M„ - J pi + 2n/£ 



cos^ enM„. - sin^ OnMnMnMn 



where 



Mn 



2n/£ Pz 



n = , 



0, n > 



yjpl + 2n/P \ p, V2^/e^ 
The solution of the above equations is given for any value of n by 



(26) 



(27) 



tan (fir, 



, tan ZOj-, 



m 



(28) 
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The vacuum state in the magnetic field reads in this case 



i^^cfl 16). 



n,p 



6t - 



(29) 



Let us now evaluate the fermion condensate b(O|'0'0|O)b. In the limit when m 
goes to zero we obtain 

lim BWmB = - lim / ^ (cos 2^0 + 2 £ cos 2^ J 



m 



+ 2E 



\eB\ 
— hm — — m 
m-^O 47r 



^ , (30) 



In — + C>(m° 



where an ultraviolet cutoff A has been introduced in order to regularize the 
integral. In contrast to the (2+l)-dimensional case, where a nonvanishing con- 
densate appears in the limit when m — > (see Eq. (22)), it is seen from 
Eq. (30) that this quantity vanishes in the (3+l)-dimensional case. Neverthe- 
less, Eq. (30) indicates that a condensate in 3+1 dimensions could in principle 
be generated in the presence of a small dynamical mass. 



4 A test example: the NJL local kernel 



At this stage it is convenient to test our formalism against the well-known 
results of the NJL model. This test, which is nontrivial, will then act as a 
validation criteria for the present formalism. 

The Lagrangian density for the original NJL model in an external field reads 



C = i){x) [i-i>'D^ - m\ tl^ix) + G {i){x)'il){x)f + {i){x)i-f^'il;{x)f 



(31) 



Let us start by discussing the (2-|-l)-dimensional NJL model. We recall that 
in 2-1-1 dimensions the 75 matrix in the chiral representation is defined as 
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where 1 is the 2x2 unit matrix. The contribution of Hint to : H2 : is easily 
obtained from Eq. (12), 



[^(1) + ^(2) _ (5(1) _ 5(2))] ^^ + (1^2), 
where the fermion propagators are given by (cf. Eq. (14)) 



(33) 



^« = — y 

iL ^ 



(34) 



Following the procedure of section 2, the diagonalization of the Hamiltonian 
: H2 : through Valatin-Bogoliubov transformations yields the following mass 
gap equations: 



(m + p) sin 200 = , n — , 

(m + p) sin 26'n - V^n cos 26'„ = , n > , 



where 



G_ 



008200 + 2 J2cos29n] . 



n=l 



The solution of these equations is 



tan 20n — 



2n 



[m + p)l ' 

which when substituted back into Eq. (36) yields 

m + p 



^1 ^2n/B + (m + p)2 ' 



Using the identity [14] 



00 -1 
2E 

n= 



/3- 



-(2n + /3)^ r(//)yo 
we can recast Eq. (38) in the form 

7r^/^£ , ^ f"^ ds 



dssi'-^e-^' coths 



G 



'-p=(m + p) r 4^e-(™+'')'^'^ coth s , 
Jo s^l^ 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



which reproduces the result obtained in leading order of the 1/N expansion [2]. 
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In the weak coupling regime, i.e. when G i , one expects the lowest Landau 
level to dominate. In this case, Eq. (38) implies 



TT 



Recalling that {0\ipip\0) = -p/{2G), we find (0|V^^|0) = -|e5|/(27r). We 
notice that this result coincides with the value of the condensate previously 
obtained in Eq. (22), in the absence of the NJL interaction. 

Let us now consider the (3+l)-dimensional NJL model. In this case, the con- 
tribution to :if2: is calculated from Eq. (13) and one obtains 

V'^ [To'S'to - 7oTr (5'7o) - ToTs-S'toTs + ToTsTr (5'7o75)] V' , (42) 
where, according to Eq. (15), the fermion propagator is 

^ = E ttV E h"(0«L(0 - vsniOvliO] . (43) 



Once again, : H2 : is diagonahzed by successive Valatin-Bogohubov trans- 
formations. For the sake of comparison with the results obtained using the 
1 /A^-expansion technique, we shall assume that pions do not condensate. In 
this case, we end up with the following mass gap equations: 



(m -|- p) sin 26*0 — {pz + 1^) cos 26*0 = , n = , 

(m + p) sin 2^„M„ - ^pl + 2n/P (cos^ ^„M„ - sin^ Or^M^MnM^) 

COS^ 9n(Tl + Sin^ ^n(sin 2(/7„(73 — COS 2</7„(7i) = , n > , 



—K 



where Mn and M„ are defined in Eqs. (25) and (27), respectively; 



^ ^ J ~2n \^ 2610 + 2 ^ sin 26'„ cos 



(44) 



(45) 



The solution of these equations is given by 



'""^" = (^:t4^' '""'^" = — ^p — ■ ^''^ 



11 



Thus, from Eq. (45) we find 



27r2£2 r m + p m + p , 



wfiich can be recast in tfie integral form 

= (m + p) r ^e-^rn+p)H^s ^^^^ ^ ^ (48) 
G Jo s 

througli the identity (39). 

In the weak couphng regime, when the LLL dominance approximation remains 
vahd, from Eq. (47) we obtain in the hmit m ^ 0, 

-^p = pln — , (49) 

where A is an ultraviolet cutoff which regularizes the integral in the momentum 
Pz- This result is consistent with the one previously found in the absence of 
the NJL interaction term [cf. Eq. (30)] and agrees with the results obtained by 
the use of other approaches [2]. Wc also notice that Eq. (49) gives the correct 
behavior near the singularity at G = 0, 

On the other hand, as G increases the expression (50) gets less accurate. This is 
due to the fact that for larger values of G all the Landau levels give important 
contributions to the dynamical mass and the LLL approximation is no longer 
valid. 

It is worth emphasizing that in the local NJL model the nontrivial chiral 
condensate is generated in the presence of an arbitrary small external magnetic 
field for arbitrary small values of the coupling constant. It turns out, as we 
shall see shortly, that this picture ceases to be valid when, in the presence of 
a homogeneous magnetic field S, we go from a local to a nonlocal fermionic 
kernel. In this case we cease to have a weak coupling limit no matter how 
small we set the kernel strength. 



5 Nonlocal kernels: the hcirmonic oscillator case 



The (2+l)-dimensional case of a nonlocal quartic kernel (given in this example, 
for the sake of simplicity, by a harmonic-oscillator kernel) constitutes, in the 
presence of a magnetic field, the first nontrivial application of our formalism. 
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The harmonic-oscillator quartic interaction is given by 



d^xd^y^P{x)-io—^p{x) -^K\x -y\^ i){y)-^^—^p{y) , (51) 




int 



where K is the interaction constant and A" are the Gell-Mann matrices, which 
account for the color structure of the interaction. This structure prevents 
the appearance of tadpole-like terms. Notice that Hint has a central (polar) 
symmetry. 

In the absence of a magnetic field, the harmonic-oscillator nonlocal NJL Hamil- 
tonian can be readily diagonalized through a "2D polar-symmetric" Valatin- 
Bogoliubov canonical transformation. The other limit, i.e. having a magnetic 
field and no quartic fermion interaction, has been already studied in the 
present paper and leads to a condensate with nonpolar, axial-like symme- 
try. It is therefore clear that we cannot accommodate both symmetries si- 
multaneously: if we had started from a sufficiently strong magnetic field, and 
progressively turned on a quartic nonlocal interaction, there will be a value 
of the quartic kernel strength where the fermion condensate cannot hold any 
longer and the system becomes disordered. 

The contribution to the Hamiltonian term :H2: in Eq. (12) is represented by 
the diagrams A and C of Fig. 1, 



Here, S'^^^x^y) is the propagator defined in Eq. (14). The tadpole diagrams B 
and D do not contribute since the A" matrices are traceless. Moreover, there 
is neither a contribution from diagram C to the '^i term nor a contribution 
from diagram A to the -02 term, because the 70 vertex is block diagonal. 

Applying successive Valatin-Bogoliubov transformations and making use of 
the simple relations 



with uJn{i) defined in Eq. (6), we arrive at the following recurrence system of 
mass gap equations: 



4{x)S^'\x,y)V{\x - y\)My) + (1^2). 



(52) 



(e - d^) iJniO = V2(n+l)u;„+i(0 , 

(e + a^)u;„(0 = V^^n-i(0, 



(53) 




(cos 26*0 + cos 2^1) sin 200^0, 



n = 0, 




(54) 



+Kf [-((2n - 1) cos2^„_i + (2n + 1) cos2^„+i) sin 2^, 



+2 y/n{y/n - 1 sin 2^„_i + ^/n + l sin 2^„+i) cos 29^ =0, n > . 
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Fig. 2. Chiral angle as a function of the Landau level n for m = and Ki^ = 0.5 
in the presence of a harmonic-oscillator potential and a constant magnetic field in 
2+1 dimensions. 

We impose the boundary condition = , so that the standard expression 
for the n = spinor, MqIO = (^o(0 > 0)"^) recovered. We also note that 
Eqs. (54) enforce 6*^ — >• 7r/4 as n — >• oo . 

Using the second equation in (54) we can calculate the chiral angle 9n+i from 
6n and 9n-i- This leaves 6i as a free parameter. By varying 9i, we can find the 
solution that minimizes the vacuum energy, 

(oo \ 2 °° 

cos 2^0 + 2 ^ cos 2^„ + Kf sin^ ^0o--J2 ^ ^in 29^ . (55) 
n=l / n=l 



The nonlinear mass gap equations (54) have no obvious analytical solutions 
and wc have to solve them numerically. As an example, we present in Fig. 2 
their numerical solution for m = and K£^ = 0.5. Having obtained the chiral 
angles we can proceed to calculate physical quantities. The chiral condensate 
{ipip) was numerically shown to converge for each value of Ki^ and, therefore, 
no renormalization is needed. This is in contradistinction to the NJL local 
kernel, which needs to be regularized. The results are shown in Fig. 3. Notice 
that when K£^ ^ we recover the value of the condensate previously obtained 
without the harmonic-oscillator interaction, i.e., —27ri'^{0\ijjijj\0) = 1. When 
K£^ ~ 0.8 the chiral condensate vanishes.^ This signals the onset of a phase 

^ Taking, for instance, K = (300 MeV)^, a value of Ki^ = 0.8 would correspond to 
a magnetic field strength of B = 1.5 x 10^^ G. 
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Fig. 3. Chiral condensate as a function of the dimensionless parameter K£^ for 
m = 0, in the presence of a harmonic-oscillator potential and a constant magnetic 
field in 2+1 dimensions. 

transition due to the presence of both the magnetic field and the fermionic 
quartic kernel (in this confining harmonic-oscillator kernel). 

In the presence of a nonvanishing current fermion mass, m ^ 0, the picture 
is quite different. As an example, in Fig. 4 we present a plot of the chiral 
condensate for m = 2 MeV. Notice that the chiral condensate diverges with 
the number N of Landau levels and, consequently, a renormalization scheme 
is needed for this case. 



6 Conclusion 



In this paper we have studied the problem of spontaneous chiral symmetry 
breaking in the case of the simultaneous presence of a homogeneous magnetic 
field and a fermionic quartic interaction, for both local and nonlocal kernels 
in 2+1 and 3+1 dimensions. The operator formalism presented here is based 
on the use of a sequence of Valatin-Bogoliubov canonical transformations, 
which allowed us to construct explicitly the vacuum of the system and also to 
calculate the chiral condensate. As illustrative examples, we have considered 
the original NJL model and the harmonic-oscillator potential. In the case of 
the NJL model we were able to recover the well-known mean-field results. 
For the (2+1) harmonic oscillator, we have shown that in the presence of the 
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Fig. 4. Chiral condensate as a function of the dimensionless parameter K£^ for a 
current mass m = 2 MeV, in the presence of a harmonic-oscillator potential and a 
constant magnetic field in 2-1-1 dimensions. 

magnetic field the system (with nonlocal interaction) exhibits new properties, 
namely, the existence of a critical magnetic field Be such that for B > Be a 
chiral condensate is induced. It would interesting to extend this analysis to the 
(3-|-l)-dimensional case. However, in the latter case the mass gap equations 
turn out to be not only recurrent but also differential and, therefore, their 
solution is more involved. 
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